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Abstract. Fluids in external gravity may oscillate with frequencies characte- 
ristic of the epicyclic motions of test particles. We explicitly demonstrate that 
global oscillations of a slender, perfect fluid torus around a Kerr black hole admit 
incompressible vertical and radial epicyclic modes. Our results may be directly relevant 
to one of the most puzzling astrophysical phenomena — high (hundreds of hertz) 
frequency quasiperiodic oscillations (QPOs) detected in X-ray fluxes from several black 
hole sources. Such QPOs are pairs of stable frequencies in the 3/2 ratio. It seems that 
they originate a few gravitational radii away from the black hole and thus observations 
of them have the potential to become an accurate probe of super-strong gravity. 
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1. Introduction: QPOs and epicyclic frequencies 

High frequency (hundreds of hertz), quasiperiodic variations in X-ray fluxes have been 
detected in several Galactic black hole sources by X-ray telescopes on spacecraft (e.g. 
van der Klis 2000, 2005; McClintock & Remillard, 2004). They are called QPOs and 
have frequencies corresponding to orbital frequencies a few gravitational radii away 
from black holes. Of special interest are the twin peak QPOs that come in pairs of two 
frequencies in the 3/2 ratio. 

Kluzniak, Abramowicz (2000, 2001) and their collaborators developed a model of 
twin peak QPOs in terms of a non-linear resonance in epicyclic oscillations of nearly 
geodesic motion of matter close to the black hole (see Abramowicz 2005 for information 
and references). In this paper we extend the model by studying epicyclic oscillations 
of relativistic fluid tori. We prove that in an axially symmetric stationary spacetime, 
an equilibrium perfect fluid torus with small cross-section always admits global epicyclic 
eigenmodes of oscillations. In Newtonian theory, the same was proved recently by Blaes 
& al. (2005). We follow closely their method, adapting it to general relativity. We 
start with a short reminder of how to construct an equilibrium model of the torus (after 
Abramowicz & al., 1978). 



2. Equilibrium of the torus 

The line element of a stationary, axially symmetric spacetime is given in the spherical 
coordinates (t, r, 9, <fi) by 

ds 2 = g u dt 2 + 2g t< f ) dtd(f) + g rr dr 2 + geed9 2 + g^dcf 2 . (1) 

We use the — h ++ signature. The metric coefficients depend on r and 9, but not on 
t and 0, and for this reason it is convenient to introduce two Killing vectors r] a = 5 a t 
and £ a = 5^. Note that, obviously, rfr} a = g tt , £ a r] a = g t4> , £ Q £ Q = g H . We assume in 
addition, reflection symmetry, i.e. that all odd ^-derivatives of the metric coefficients 
vanish in the symmetry plane 9 = ir/2. 



2.1. Circular orbits and circular geodesies 

We start by briefly recalling some known results (see e.g. Wald, 1984) concerning circular 
orbits and circular geodesies. A circular orbit is given by 

u a = A(r t a + n£ a ), (2) 

with Q being the angular velocity, and 1/A 2 = —(f] a f]a + 2Qi] a C, a + f2 2 £ Q £ Q ). The 
conserved energy and angular momentum are given by £ = —r] a u a , C = (, a u a , and the 
conserved specific angular momentum by £ = £/£. The acceleration at the circular 
orbit is 

d a g tt -2M a g t <{> + e 2 d a g^ \ J_ ( &J_ \ 

g tt _ 2£gt4> + £2gM J 2U\dx a )/ U 
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where U = g u — 2£g t ^ > + £ 2 g^ is the effective potential. 

The geodesic circle is defined by the condition a a = 0. Radial u r and vertical ujq 
epicyclic frequencies for circular geodesies around black holes have been derived first by 
Kato & Fukue (1980), see also Silbergleit & al. (2001). They may be expressed in terms 
of the second derivatives of the effective potential (Wald 1984, Abramowicz & Kluzniak 
2002), 

2 _ S 2 (<nA\ 2 _ S 2 (&U\ 

U * ~ 2A 2 g rr \dr 2 ) t ' " 2A 2 g ee \d6 2 ) ' W 

These are frequencies measured by an observer at infinity. Note that these expressions 
are invariant — independent of the choice of metric signature and coordinates. 

2.2. Perfect fluid 

The torus is made of perfect fluid with stress-energy tensor T 01 ^ = (e + p)u a up + p5 a p, 
where e and p are the energy density and pressure of the flow. The energy density 
consists of rest mass density p and thermodynamic internal energy density U. 

The dynamics is governed by the conservation law V^T^ = 0. Projection of this 
into the 3-space perpendicular to the 4- velocity, together with the 4-acceleration formula 
©, gives 

- ^ - k 2 v„« = 

e + p 2 

l -(8 2 - 8l)W a g tt - (£ 2 £ - ^ )V a ^ + \{E 2 £ 2 - «)V Q ^. (5) 

From (J3J) it follows that the pressure maximum V Q p = is a geodesic circle a a = 
and lies in the equatorial plane cos 8 = 0. We shall call this circle the torus center and 
indicate by subscript quantities evaluated there; in particular £ , £ are the energy 
and specific angular momentum at the torus center. 

2.3. Barotropic Fluid 

For a barotropic fluid p = p(e), the left hand side of this equation is a gradient of a scalar 
function. Thus, the right hand side must also be a gradient of a scalar function. We 
denote it by V a $. Obviously, both sides of equation © are gradients in the particular 
case of a relativistic polytrope (which we now assume), i.e. when e = np + p and 
p = Kp 1+l / n . In this case the enthalpy H equals, 

dp 



H 



p 

; — — - = In (1 + n) 

' p 



l + n)- + l 



1 + (6) 



[1 + n)p + (p/K) n l n+1 

The last approximation is valid when the local sound speed c 2 oc p/p is negligible with 
respect to the speed of light c 2 = 1. Substituting H into equation (JHJ) and integrating, 
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we find the Bernoulli equation in the form 



P 



-EqU + (1 + n)— — const 



(7) 



nc 



sO 



F 2 

2 v 



(8) 



We define the function / = f(r, 6) by, 

- = -f{r,0), f{r,0) = 
P Po 

where c 2 = (n + l)p/(np) is the square of the sound speed at the torus center. We 
now calculate the function / in the vicinity of the center of the torus. For this purpose, 
we introduce new coordinates x and y by dx = y/g~^dr/r and dy = —^/geed9/r , and 
the condition that x — = y at the torus center. The metric coefficients in the above 
definition are evaluated at the center of the torus. For the effective potential we obtain 
in terms of x, y, 

where all metric coefficients are evaluated at the torus center. In the vicinity of the 
equilibrium point the potential \I/ can be expressed as 



U 



2 



1 

Qrr 



f 2 



2g rr to 
Finally, by substitution into equation 



Or 



f=l- 



2Un 



where 



(g tt + ^o^) 2 



df 

dr 



2 2 

r n x . 



and using equation (j3J, we obtain 



2,-22 

x +uj e y 



{g tt r - 



t<t>\ 

,r J ) 



P 2 = 



(2nc s 2 ) 



(10) 



(12) 



and u r and ujq are ratios of the epicyclic frequencies to the orbital frequency f2 at the 
center of the torus. 

The torus is slender when /3 <C 1. The approximation used in (JBJ) is obviously 
correct for slender tori. It is interesting to note that slender tori with constant angular 
momentum distributions £(r, 8) = have elliptical cross-sections with semiaxes in the 
ratio of the epicyclic frequencies. The same is valid for Newtonian slender tori (Blaes & 
al., 2005). In the particular case of the Newtonian 1/r gravitational potential, 
and slender tori have exactly circular cross-section (Madej & Paczyhski, 1977). 



3. Perturbed Equilibrium and Epicyclic Modes 

3.1. Barotropic Tori and the Papaloizou-Pringle Equation 

We consider small perturbations around the stationary and axially symmetric equi- 
librium of the torus (slender or not), in the form 5X(r,9,(j),t) = 5X^{r 1 6)e 1<ym ^ U)t \ 
The equation that governs these perturbations was derived in Newtonian theory by 
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Papaloizou & Pringle (1984) in terms of a single quantity W. Its general relativistic 
version has the form 

W = » (13) 
u t p{m\l — to) 

For simplicity we assume for now that the fluid angular momentum I = const. We also 
assume that the flow is potential in the sense that 

u a = — ^V Q $, (14) 
e + p 

where $ is some scalar field. Linearizing this equation along with the normalization of 
the four-velocity, we find that 5$ = iW, so that 

S Ua = — — \iW a -u a uHmQ, -u)W] . (15) 
e + p 

Combining these with the perturbed continuity equation, V At (pw M ) = 0, 

V /t (pcW) = -i{mtt - u)u l 8p, (16) 
we arrive at the relativistic version of the Papaloizou-Pringle equation 



J d_ 

-g) 1 / 2 dx l 



2nA 2 (mVL 



- (mY* ~ 2 W + ^g tt )PW + ^ P~ l W = 0, (17) 

where g is the determinant of the metric and we used c 2 s <C c 2 . In the slender torus limit 
f3 — > 0, this equation becomes 

d 2 W dfdW d 2 W dfdW 
ox z ox ox oy z oy oy 

where x = x//3, y = y//3, and v = uj/VLq. 

Equation (|18j) is identical in form with the nonrelativistic Papaloizou-Pringle 
equation for oscillations of the constant angular momentum slender torus (cf. Blaes 
& a!., 2005), so one could take advantage of Feynman's reminder, "the same equations 
have the same solutions". In the particular case u r = ujq the nonrelativistic Papaloizou- 
Pringle equation was fully solved in exact analytic form by Blaes (1985) who gave a 
complete set of its eigenmodes, with a complete analytic description of all eigenfunctions 
and all eigenfrequencies. Blaes' solution cannot be directly applied to the relativistic 
Papaloizou-Pringle equation (fTTjl. because for relativistic tori u r ^ cug. 

More recently, Blaes & al. (2005) found that in the general Newtonian case (non 
spherically symmetric potential, u r ^ there are two particular solutions of the 
general Papaloizou-Pringle equation, 

W r = C r 5e i(m *- Wrt) W z = C z ye L{m4> ~ UBt \ (19) 

with C r , C z being two arbitrary constants. We claim that they are also solutions 
of the relativistic Papaloizou-Pringle equation 1)18)1 . One may check this by a direct 
substitution of the solution (JT§|) into equation (JTHJ). 



>2 
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The solution (jTUjl is consistent with epicyclic modes. Indeed, the fluid velocity is 
spatially constant on the torus cross-sections, entirely radial in the case of the radial 
epicyclic mode W r and entirely vertical in the case of the vertical epicyclic mode W z . 
The radial and vertical oscillations are sinusoidal, with frequencies equal to the epicyclic 
frequencies u r and ug. 



3.2. Baroclinic Tori 



We assumed previously that the tori had constant specific angular momentum and were 
isentropic. This simplifies the perturbation equations enormously by removing the local 
restoring forces that give rise to internal inertial modes. We now consider more general 
configurations, and show that global epicyclic modes are still present even here. 

The equilibrium of a relativistic baroclinic torus is described by equation (jSJ), and 
the equations for axisymmetric perturbations are those of conservation of rest mass, 
momentum, and adiabatic flow: 



t d5p tt dSu t , td,^ u 4> 
u -^7- + P9 —z h P9 



dt 



dt 



dt 



d 



-g) 1 / 2 dr 



-g) 1/2 p5u r 



u 



t d5ui 



Ot 



£ 



9a 



9rr 

o t( i> 



d 



(- g y/ 2 de _ 



-g) 1/2 p6u e 



9ee 



(20) 



+ 



u l g 



P,i 



e + p) 2 
dt 



{6e + Sp) 



Sut + { !h 
1 



t<t> 



Igf) 5u 4 



e + p 



(21) 



d5u t 



and 



^3 [g^g^-^f] 5u-VZ 

e + p\ d5p 



+ gW ddp 
e + p dt 



tdSp t 2 



P 



dt 



+ Su ■ Vp — c 



e + p 



5u ■ Vp = 0. 



Here g is the determinant of the metric, i = r or 9, and 



dp 
de 



P 



dp 



(22) 



(23) 



(24) 



e + p \dp / 

is the adiabatic sound speed. These equations form a closed system when we include 
the perturbed normalization of the four-velocity, 

Sut = —Q5u^,, (25) 

and a perturbed equation of state, 



/ de 



\dp 



Se = — 5p + — Sp. 



de\ 



d P J , 



(26) 
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The relativistic perturbation equation for rest mass conservation can be simplified by 
neglecting derivatives of the metric in the slender torus limit. Equation (J2Uj) therefore 
becomes 



u 



,dSp 



+ P9 



tdSut 



+ P9 



dt ' ™ dt ' r " dt 
We now seek solutions of these equations in the slender torus limit in which Su 
and Sug are spatially constant. Equations (J22j) . (J25j) . and (J27)l then imply that 



dSp 




dt 


M 


dS p 


i 


dt 




+ h\ 










L(u*) 



l -8u ■ Vp + (e 3 [(^) 2 - »«^] 5u ■ W 



(28) 



-Su ■ Vp 

P* 



(29) 



and 



dSu 4 

~dT 



i 95«t 

(u t v$ 



+ 9**) 



D \£u*(e+p) [(g 



t<t>\2 



Su-Vp- £ 2 Su ■ VI 



where 



D = 1 



iu 



(30) 



(31) 



Now, in the slender torus limit, c 2 , pj p, and pj e are all very small. Gradients in pressure 
and density also scale as one over the small thickness of the torus. Hence D — > 1 and 
equations (J28|) - (j30|) reduce to 
d5p I 



and 



dt 

dSp 
~dt 

dSu d 



;SU ■ Vp, 



— -Su ■ Vp, 

u 



1 85ut 



-S 2 5u ■ V£. 



dt Q dt 

Differentiating equation (J2T|) with respect to time, and using equations 
(jH2J)-((HlJ) to eliminate Se, Sp, Sp, Su t , and Su^,, we obtain 



(32) 
(33) 

(34) 
and 



dt 2 



P.i 



(e + p) 



1 du l Vp 



u l dx % 



e + p 



+ u t £ 3 (Su-V£) 



t<f>\ 



9,1 + 



(35) 
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The second term on the right hand side of this equation is negligible compared to the 
others in the slender torus limit. Using the fact that the four-acceleration vanishes at 
the pressure maximum, the following relationships are true in the slender torus limit: 



u t £ 3 QV£ = ]-V£ 2 = \v{£ 2 - £ t 



2 ) 



u t £ i £V£ 



1 



2n2\ 



-V(£ 2 £ 2 -£ 2 > 



and 



m'£ 3 (1 + ttl)\7£ = V{£ 2 £) = V{£ 2 £ - £ 2 £ ) 
Hence equation (|35|) may be rewritten as 



(36) 
(37) 

(38) 



u 



t\2 



d 2 5ui 
dt 2 



Su-V 



P.i 



e + p) 
u-V(£ 2 



+ ±g*6u ■ V(£ 2 - £ { 



+ -gfSvi ■ V(W - S, 



" 2 f 2 ) 



(39) 



In the slender torus limit, the metric derivatives may be absorbed inside the gradients, 
so that 

,d 2 5iij 



u 



1 



P.. 



1 



L(e + ,) + 2^ 2 ^ 



gf(£ 2 £-£Z£o) + ^9T(£ 2 t 2 -tftl) 



1 



£$6u- V(VW)| { 



(40) 



where the last equality comes from the equilibrium equation ([5]). Comparing equation 
PUj) with equation (jlj), we see that we once again have radial and vertical modes which 
oscillate at the respective relativistic epicyclic frequencies. 



4. Discussion and conclusions 

We have shown that the equations of motion of a fluid slender torus in axially symmetric, 
stationary spacetimes admit small oscillations occuring at the epicyclic frequencies |. 
For an initially axisymmetric torus, these two oscillations are exact solutions of the fluid 
perturbation equations in the limit of vanishing sound speed (when compared to the 
orbital frequency), and therefore represent two modes of oscillation of a slender torus. 

| This is a specific example of a more general property. Consider an isolated fluid body, moving in a 
fixed spacetime. Synge (1960) proved that inside the world tube of the body a geodesic line exists. For 
the fluid tori considered here, this is the geodesic circle at the pressure maximum. A small perturbation 
of the body results in a perturbed world tube that also contains a geodesic line, close to the original one. 
The spacetime distance between them is given by the geodesic deviation equation that in the special 
case of a geodesic circle © in the spacetime {Q reduces to two uncoupled epicyclic oscillations with 
frequencies |@J|. Global oscillations of fluids can have frequencies slightly different from the geodesic 
epicyclic frequencies, as equation (|41|l shows for nonslender tori. See also Biesiada 2003, Colistete, 
Leygnac & Kerner 2002, and Kerner, van Holten& Colistete, 2001. 
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One of these modes corresponds to a uniform vertical displacement of the torus at the 
vertical epicyclic frequency, the other to a purely radial motion at the radial epicyclic 
frequency. 

In this paper (Paper II) we considered linear deviations from circular geodesic 
motion. In the linear approximation the two epicyclic modes are uncoupled and 
independent. In reality, non-linear terms would inevitably couple the two modes. In 
the other papers of this series we discuss how the coupling leads naturally to the 3/2 
resonance. Because the 3/2 resonance is the main motivation for the whole series, we 
summarize here an idea as to how the resonance is excited and why the ratio 3/2 is a 
natural consequence of strong gravity. 

For Newtonian slender tori, Paper III shows that (slightly) away from the exact 
slender torus limit (3 = (which corresponds to c s — 0), the frequencies of vertical and 
radial oscillations of fluid are (see also Lee & al., 2004), 

(u$ f = u 2 - A z (n) c 2 s , K 2 ) f = uj 2 - A r (n) c 2 s , (41) 

where A z (n) and A r (n) are explicit positive functions of the polytropic constant n. In 
Paper I (Kluzniak & Abramowicz, 2002), the same result was derived from very general 
arguments, but without the specific form of A z (n) and A r {n). In Paper I we also argued 
that the presence of c 2 in both frequencies provides a pressure coupling between these 
two modes and that, accordingly, one may consider the vertical mode as being described 
by the Mathieu type equation, 




Figure 1. Resonance conditions in Schwarzschild spacetime (from Torok & al., 2005) 
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It is well known that a parametric resonance occurs when 

^ _ 2o^ n— 1,2,3,... (43) 
n 

For black holes, ujq > u r , and therefore the smallest value of n consistent with the 
resonance is n = 3 (Paper I), which corresponds to the 3/2 ratio of the vertical and 
radial epicyclic frequencies. This is shown in Figure ^ along with a few other possible 
rational ratios of the vertical and radial epicyclic frequencies. 

The radius (and frequency) at which the 3/2 resonance occurs depends on the black 
hole spin a. Therefore, for sources with known mass, the observed QPO frequencies may 
be used for a direct measurement of the black hole spin, as was first done by Abramowicz 
& Kluzniak (2001). 
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